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Chandra observation of two shock fronts in the merging galaxy 
cluster Abell 2146 
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ABSTRACT 

We present a new Chandra observation of the galaxy cluster Abell 2146 which has revealed a 
co mplex merging system w ith a gas structure that is remarkably similar to the Bullet cluster 
(eg. iMarkevitch et al.ll2002h . The X -ray image and temperature map show a cool 2 — 3 keV 
subcluster with a ram pressure stripped tail of gas just exiting the disrupted 6 — 7 keV primary 
cluster. From the sharp jump in the temperature and density of the gas, we determine that the 
subcluster is preceded by a bow shock with a Mach number M = 2.2 ± 0.8, corresponding to 



a velocity v = 22QiQi^'^^ km s^' relative to the main cluster We estimate that the subcluster 
passed through the primary core only 0.1 —0.3 Gyr ago. In addition, we observe a slower 
upstream shock propagating through the outer region of the primary cluster and calculate a 
Mach number M = 1 .7 ± 0.3. Based on the measured shock Mach numbers M ^ 2 and the 
strength of the upstream shock, we argue that the mass ratio between the two merging clusters 
is between 3 and 4 to one. By comparing the Chandra observation with an archival HST 
observation, we find that a group of galaxies is located in front of the X-ray subcluster core 
but the brightest cluster galaxy is located immediately behind the X-ray peak. 
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1 INTRODUCTION 

Galaxy clusters are assembled by hierarchical mergers of smaller 
subclusters and groups. These subclusters collide at velocities of 
a few thousand kms^', releasing as much as lO^'* erg of ki- 
netic energy as thermal energy by driving shocks, generating tur- 
bulence and probably accelera ting relativistic particles (see eg. 
ISarazinll200ll : lFeretti et al.ll2002ll . Major cluster mergers are there- 
fore the most energetic events since the Big Bang. Shock fronts 
provide a key observational tool in the study of these systems. 
They can be used to determine the velocity and kinematics of the 
merger and to study the conditions and transport processes in the 
IC M, including electr on-ion equilibrium and thermal conduction 
(eg. lMarkevitchll2006h . Combining X-ray observations of merging 
clusters with gravitational l ensing st udies has also produced direct 
detec tions of dark matter Jciowe et al . 2004, 2006; Bradac et al. 
l2006h and constraints on the dark matter self-interaction cross 
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section jMarkevitch et alj |2004| : iRandall et all |2008|) . Radio ob- 
servations of supernova remnants indicate that a fraction of the 
shock energy can be converted into the acceleration of relativis- 
tic particles (eg. Blandford & Eichler 1987). It is likely that this 
process also operates in cl uster mergers and could produce syn- 
chrotron radio emission (e g.lFeretti et alj2002l : lFerettil2005l : iBuotd 
I2OOII : lKempner&Davidll2004h and inverse Compton hard X- 
ray em ission (eg. Fusco-Femian o et"al]| 19991 : iRephaeli et al.|[l999l : 
iFusco-F emiano et al. 2005). However, as radio emitting electrons 
have short radiative lifetimes (10^ — 10** yr), it is difficult to ex- 
plain the ~ M pc size of extended radio halos (for a review see eg. 
lBrunettill2003l) . 

X-ray observations of merging shocks currently provide the 
only method for deter mining the velocity of t he cluster gas in the 
plane of the sky (eg. IMarkevitch et al.lll999l) . By measuring the 
temperature and density of the gas on either side of the shock 
using X-ray imaging spectroscopy, the shock velocity can be cal- 
culated from the Rankine-Hugoniot jump conditio ns. While many 
clusters are found to have shock-heated regions (eg. 
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the detection of a sharp density edge and an unambiguous jump 
in temperature is rare. Currently only two shock fronts have 
been found by Chand r a: one in the Bul let cluster (IE 0657-56; 
iMark evitch et al."2002': Marke vitchlllOOel) and the other in Abell 
520 (^Markevitch et al. 2005). In this paper we present two new 
merger shock fronts discovered in a recent Chandra observa- 
tion of the gal axy clus ter Abell 2146 at a redshift z = 0.234 
dStruble & Roo d 1999; B ohringer et"ai]|2000h . 

We assume //o = 70 km s"' Mpc"' , a„, =0.3 and Sl^ = 0.7, 
translating to a scale of 3.7 kpc per arcsec at the redshift z = 0.234 
of Abell 2146. All errors are IcT unless otherwise noted. 



2 DATA PREPARATION 

Abell 2146 was observed with the Chandra ACIS-S detector for 
43 ks split into two observations which were taken only a day apart 
in April 2009 (Obs. IDs 10464 and 10888). The data were anal- 
ysed with CIAO version 4.2 and CALDB version 4.2.0 provided 
by the Chandra X-ray Center (CXC). The level 1 event files were 
reprocessed to apply the latest gain and charge transfer inefficiency 
correction and then filtered to remove photons detected with bad 
grades. The improved background screening provided by VFAINT 
mode was also applied. The background light curve extracted from 
the ACIS-S 1 level 2 event file was filtered using the LC.CLEAN 
scripQ provided by M. Markevitch to identify periods affected by 
flares. There were no flares in either observation of Abell 2146 so 
we proceeded with the final cleaned exposure of 43 ks. 

As the two separate observations were taken so closely to- 
gether, with effectively identical chips positions and roll angles, 
we were able to reproject them to a common position (Obs. ID 
10464) and combine them. Exposure-corrected images were cre- 
ated by combining the two cleaned event files and assuming a mo- 
noenergetic distribution of source photons of 1.5 keV, which is ap- 
proximately the peak energy of the source. 



3 IMAGING ANALYSIS 

Exposure-corrected images of the galaxy cluster and a zoom in of 
the cluster core are shown in Fig. [T] The X-ray emission is ex- 
tended from SE to NW and appears to cut off abruptly at either 
end of this axis. The bright, dense core has been displaced from 
the cluster centre and is being stripped of its material to form a tail 
of gas towards the NW. The brightest cluster galaxy (BCG) shown 
in Fig. [2]is displaced ~ 10 arcsec to the NW of the X-ray surface 
brightness peak and contains a point source detected in the Chan- 
dra X-ray image and a VLA 4.9 GHz radio image (NRAOA'LA 
Archive Survey). The X-ray point source is detected at hard ener- 
gies (2—10 keV) and it is likely that this corresponds to an AGN, 
but a measurement of the flux was difficult with the superimposed 
cluster emission. We estimated the point source flux by extracting 
the source counts in a region of 2 arcsec radius and subtracting 
the cluster emission using a surrounding region from 2 — 5 arcsec 
radius. Using a powerlaw model with photon index F = 2, we esti- 
mated the point source luminosity in the energy range 2—10 keV 
to beL2_iokeV = 1.6 ±0.4 X 10"*^ erg s"' . 

The unsharp-masked image shown in Fig. [T] (right) highlights 
the different structures in the cluster gas. The cluster is separated 
into two concentrations of X-ray emitting gas: the cluster core and 
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its trail of emission to the SE and a second, more diffuse region 
of material to the NW. Figs. [T] and [3] also reveal several edges in 
the X-ray surface brightness. The sharpest of these defines the SE 
edge of the bright subcluster core. A second outer edge is visible 
~ 0.5' to the SE of the core although the steep decline in the cluster 
surface brightness reduces the significance of this feature. The third 
edge appears at the NW edge of the galaxy cluster. There could 
potentially be a fourth surface brightness edge separating the two 
concentrations of cluster gas. Finally, the subcluster's tail of gas 
also appears to have an extension to the SW. 

Abell 2146 a ppears to have a remarkably similar struct ure to 
the Bullet cluster i lMarkevitch et alJl2002l ; lMarkevitchll2006l ). The 
X-ray morphology suggests a recent merger where a subcluster 
containing the dense core has passed through the centre of a second 
cluster, the remnant of which appears as the concentration of gas to 
the NW. The dense subcluster has just emerged from the primary 
core, travelling to the SE, and is trailing material that has been ram 
pressure stripped in the gravitational potential. The passage of this 
dense subcluster core is likely to generate sharp shock fronts in the 
ICM. 

To confirm the detection of the edges in the X-ray emission, 
we produced surface brightness profiles in two sectors to the NW 
and SE from the AGN, as shown in Fig.|4] The surface brightness 
profiles were extracted from the merged exposure-corrected image 
in the energy range 0.3 — 7.0 keV. The radial bins are 1 arcsec wide 
in the cluster centre and then increase in size as the cluster surface 
brightness declines and the background subtraction becomes more 
important. Point sources were identified using the CIAO algorithm 
WAVDETECT, visually confirmed and excluded from the analysis 
using elliptical apertures where the radii w ere conservatively se t 
to five times the measured width of the PSF jFreeman et a i]|2002h . 
The background was determined in a sector taken from the SE edge 
of the ACIS-S3 chip, 200 - 270 arcsec from the central AGN, in a 
region that is largely free of cluster emission (Fig.|4}. 

Fig. [5] (left) shows the SE surface brightness profile, centred 
on the AGN. There is steep decline in the surface brightness at 
18 arcsec radius, marking the edge of the subcluster core, and then 
a second break at 55 arcsec. Fig. [5] (right) shows the NW surface 
brightness profile declines slowly through the length of the ram 
pressure stripped tail to a break at ~ 50 arcsec, marking the sep- 
aration between the two concentrations of X-ray gas. There is a 
sharp edge in the surface brightness profile at ~ 120 arcsec, where 
the surface brightness drops by a factor of ~ 7. In summary, we 
identify four surface brightness edges: two in the SE sector at radii 
of 18 arcsec and 55 arcsec in front of the subcluster core, and an- 
other two in the NW sector at radii of 50 arcsec and 120 arcsec. 
To determine whet her these surface brightness edges are shocks 
or cold fronts (eg. IMarkevitch et al.l l 2OO0l ; Ivikhlinin et ai]|200ll ; 
IMarkevitch & Vikhlininll2007h , we extracted and analysed X-ray 
spectra on either side of these features to determine the ICM tem- 
perature and density. 



4 SPECTRAL ANALYSIS & DEPROJECTION 

We first extracted an overall cluster spectrum using an outer ra- 
dius of 2', which contained the vast majority of the cluster emis- 
sion, after excluding point sources. The background was sub- 
tracted using a spectrum extracted from a cluster-free region at 
the edge of the chip and appropriate responses and ancillary re- 
sponses were generated. The spectrum was restricted to the en- 
ergy range 0.5 — 7.0 keV and grouped to contain a minimum 
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Figure 1. Left: Exposure-corrected image in the 0.3-5.0 keV energy band smootlied with a 2D Gaussian fj = 1.5 ai'csec (North is up and East is to the left). 
The logarithmic scale bar has units photons cm^- s^' arcsec"- . Right: Unsharp-masked image created by subtracting images smoothed by 2D Gaussians 
with a = 5 and 20 arcsec and dividing by the sum of the two images. 




Figure 2. Hubble Legacy Ar chive image of the brightest cluster galaxy in " 3 " 2 -o 1 o o.i 0.2 3 

Abell 2146 jSand et 3112001 with Chandra X-ray and VLA 4.9 GHz radio 

(NRAOA^LA Archive Survey) contours superimposed in red solid and blue Figure 3. Unsharp-masked image as in Fig.[T]with structural features in the 

dashed lines respectively. cluster labelled. 



of 30 counts per spectral channel. A single temperature fit to 
the cluster spectrum usi ng an absorbed thermal plasma emission 
model PH ABS(MEKAL) (Balucinska-Church & McCammon' 19921; 
iMewe et al. J985; 1986; Kaastra 1992; Liedahl et al. 1995) in 
XSPEC version 12.5.0 dAmaudl Il996h produced a good fit (re- 
duced = 0.99 for 328 degrees of freedom) with a temper- 
ature keV, luminosity Lx = 1.55 ±0.02 x lO'*^ ergs"' 



(0.01 - 50.0 keV) and an abundai ice of 0.37 ± 0.04 Z,^ mea sured 
assuming the abundance ratios of lAnders & Grevessd ( Il989h . The 
Galactic hydrogen column density was left as a free parameter, 
giving a value rin = 0.028 ±0.004 x 10^ ^ cm~^ which is con sis- 
tent with the Galactic value measured bv lKalberla et alj ( l2005h of 
nu = 0.03 X 10^2 cm-2 . 

We also fitted an absorbed single temperature model to a spec- 
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Figure 4. Exposure-corrected image of the ACIS-S3 chip (8x8 arcmin) in 
the 0.3-7.0 keV energy band smoothed with a 2D Gaussian a = 2.5 arcsec. 
The sectors and background region used to produce surface brightness pro- 
files are labelled. The dashed lines across the sectors mark the approximate 
location of the surface brightness edges. 

trum extracted in a region of radius 10 arcsec, which approxi- 
mately encloses the bright subcluster core, and excluded the AGN. 
The hydrogen column density parameter was poorly constrained 
and so was fixed to the Galactic value. This produced a temper- 
ature of 2.7 ±0.1 keV and a metallicity of 0.6±0. IZ© (reduced 
= 1.4 for 68 de grees of freedom). The spectral fit was signif- 
icantly improved by adding a cooling flow MKCFLOW component 
which models gas cooling down to low temperatures {x^ = 94 re- 
duced to = 80 for 67 degrees of freedom). The low tempera- 
ture limit of the MKCFLOW model was fixed to 0.1 keV and the 
higher temperature and metallicity were tied to the MEKAL com- 
ponent parameters. The MKCFLOW model normalization suggests 
that 40 ± IOMq yr^^ could be cooling out of the X-ray and down 
to low temperatures in the subcluster core. 

If we exclude the cool core from the overall cluster spec- 
trum we get a cluster temperature of 7.5 ±0.3 keV and luminos- 
ity Lx = 1.55 ±0.02 X 10"*^ erg s"' (0.01 - 50.0 keV). Therefore 
Abell 2146 falls on the L^ — T relation for local clusters jPratt et aEI 
l2009l) . 

4.1 Temperature and density maps 

We used spatially resolved spectroscopy techniques to produce 
maps of the projected gas properties in the cluster core (Fig. [6](. 
The central ~ 4 x 4 arcmin was div ided into bins using the Contour 
Binning algorithm ( ISandersll200^ . which follows surface bright- 
ness variations. Regions with a signal-to-noise ratio of 32 (~ 1000 
counts) were chosen, with the restriction that the length of the bins 
was at most two and a half times their width. The background spec- 
trum was subtracted from the observed dataset and appropriate re- 
sponses and ancillary responses were generated. The spectra were 
grouped to contain a minimum of 20 counts per spectral channel 
and restricted to the energy range 0.5 — 7.0 keV. Each spectrum was 
fitted in XSPEC with an absorbed MEKAL model with the absorption 
fixed to the Galactic value iin = 3.0 x 10^" cm^^ jKalberla et al.l 



I2OO5I) and the redshift fixed to 0.234. The fitting procedure min- 
imised the ;|j^-statistic. The errors were approximately ~ 15% in 
temperature and '-^ 8% for the emission measure. However, the high 
temperature bins greater than 10 keV are poorly constrained by 
the energy range of Chandra producing larger errors greater than 
~ 30%. We fixed the metallicity to an average value of 0.4Zq in 
the spectral fits; the limited number of counts produced a poorly 
constrained metallicity parameter if it was left free. However, we 
found that fixing the metallicity still produced a very similar tem- 
perature and density map. 

Shown in Fig. [6] are the projected emission measure per unit 
area, temperature and 'pressure' maps. The projected 'pressure' 
map was produced by multiplying the square root of the emis- 
sion measure per unit area and the temperature maps. The pro- 
jected emission measure map shows the strongly peaked core sur- 
face brightness, bright tail of stripped material and the elongated 
morphology of the cluster in the NW to SE direction. There could 
also be a spur of emission out to the SW of the subcluster core. The 
projected temperature map shows strong variations across the clus- 
ter. In the dense subcluster core, the temperature drops as low as 
1 .9 ± 0. 1 keV and then sharply increases to the SE up to 8 — 10 keV 
between the two SE surface brightness edges. The large errors on 
the high temperature 10^2 keV bin immediately in front of the cool 
subcluster make it consistent with a constant temperature between 
the edges. The ram pressure stripped tail of material appears as a 
warmer stream of gas (5 — 8 keV) behind the subcluster core and 
trails back to the hottest region of the disrupted main cluster. The 
SW spur of emission suggested in the emission measure map cor- 
responds to a 5 — 7 keV region of gas which could be connected to 
the ram pressure stripped tail from the core. The high temperatures, 
lll^keV, intheNW region of the cluster suggest shock heating of 
the gas is likely. The sudden drop in temperature at the NW edge 
of the cluster, which coincides with the surface brightness edge, is 
particularly suggestive of a shock front. 

The NW edge is clearly visible in the projected 'pressure' 
map (Fig. [6] right) as a sudden increase and subsequent drop in 
the pressure by a factor of ~ 5. To the SE of the subcluster core, 
the pressure is approximately constant for ~ 20 arcsec and then 
drops abruptly at larger radii. This could indicate that a contact 
discontinuity or cold front immediately precedes the subcluster 
rather t han a shock ([M arkevitch et al...2000. : .Vikhlinin et al...200ll: 
iMarkevitch & Vikhlininll2007h . 

4.2 Projected & deprojected radial profiles 

We obtained a more significant detection of the temperature and 
density changes across the surface brightness edges by extracting 
radial profiles in the NW and SE sectors (Fig.|7j. The radial bins 
were positioned so as to determine the gas properties on either side 
of the surface brightness edges while maintaining a minimum of 
3000 counts in each extracted spectrum. This lower limit ensured 
enough counts to provide a good spectral fit and constraints on 
the parameters. The steep decline in the cluster surface brightness 
in the outermost radial bins made it particularly difficult to con- 
strain the outer gas properties. A wide outer radial bin from 1 to 
2.5 arcmin was therefore required for the SE sector to determine 
the nature of the surface brightness edge in front of the dense core. 

Spectra were extracted from the radial bins in each sector. 
Point sources were identified and excluded as detailed in section 
[3] The spectra were analysed in the energy range 0.5 — 7.0 keV 
and grouped with a minimum of 30 counts per spectral bin. The 
background subtraction was particularly important for the regions 
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Figure 5. Surface brightness profiles in the energy range 0.3-7.0 keV for the SE (left) and NW (right) sectors as shown in Fig.|4] The approximate locations 
of the surface brightness edges shown in Fig.|4]are marked with dashed Unes. Both profiles are centred on the AGN. 
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Figure 6. Left: projected emission measure per unit area map (units cm" arcsec" ). The emission measure is the XSPEC normalization of the MEKAL 
spectrum K = EI/ (4 x 10'*;rD^(l +z)^), where EI is the einission integral EI = J n^nudV. Centre: projected temperature map ( keV). Right: projected 
'pressure' map (units keV cm"' arcsec"^ ) produced by multiplying the emission measure and temperature maps. The blue lines correspond to the dashed 
lines in Fig. [5] The excluded point sources are visible as small white circles. 



outside the surface brightness edges where the flux is low. We com- 
pared the use of an on-chip background extracted from a source- 
free sector at the edge of the chip (Fig. with a blank-sky back- 
ground extracted from the data sets available from the CXC and 
normalized to the count rate in the high energy band 9.5 — 12 keV. 
Although the use of a blank-sky background would account for 
spatial variations in the background count rate across the chip, 
this was estimated to be only a few per cent in the energy band 
0.5 — 7.0 keV. The differences between the blank-sky and on-chip 
backgrounds when extracted from the same source-free region were 
found to be more significant. The normalized blank-sky spectrum 
overestimated the background at low energies < 1 keV and under- 
estimated the background count rate in the energy range 2 — 5 keV 



by a factor of ~ 1.3 which significantly affected the temperature 
result in regions of low surface brightness. We therefore used a 
background spectrum extracted from the cluster-free sector at the 
edge of the chip for the spectral analysis. Response and ancillary re- 
sponse files were generated for each spectrum, weighted according 
to the number of counts between 0.5 and 7.0 keV. These projected 
spectra were then fitted in XSPEC with an absorbed single tempera- 
ture thermal plasma emission model PHABS(MEKAL). The redshift 
was fixed to z = 0.234 and the absorption was fixed to the Galactic 
value riH = 0.03 x 10^^ cm"^ jKalberIa et ai]|2005h . 

However, these projected spectra extracted from the cluster on 
the plane of the sky correspond to summed cross-sections of all the 
emission along the line of sight. A spectrum from the centre of a re- 
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Figure 7. Exposure-corrected image in the 0.3-7.0 keV energy band 
smoothed with a 2D Gaussian a = 2.5 arcsec. X-ray surface brightness con- 
tours have been overlaid (thin red solid hues). The regions used to analyse 
the NW and SE surface brightness edges are shown by thick black solid and 
dashed lines respectively. 

laxed cluster will therefore contain a range of spectral components 
from the core to the cluster outskirts. To determine the properties of 
the cluster core these projected contributions from the outer clus- 
ter layers should be subtracted off the inner spectra by deprojecting 
the emission. However, deprojection routines require information 
about the line of sight extent of the cluster and generally assume 
that the cluster is spherically symmetric. While this is a reasonable 
assumption for a relaxed cluster, the highly irregular and elongated 
morphology of Abell 2146 clearly deviates from spherical sym- 
metry. We have therefore compared the projected and deprojected 
spectra and discussed the validity of the assumption of spherical 
symmetry for each sector considered. 

We used a straightforw ard model-independent spectral depro- 
jectio n routine (DSDEPROJ: ISanders & Fabianir2007l : iRussell et al.1 
which assumes spherical symmetry. DSDEPROJ starts from 
the background-subtracted spectrum extracted from the outermost 
annulus and assumes it was emitted from part of a spherical shell. 
This spectrum is scaled by the volume that is projected onto the 
next innermost shell (geometric factors from^ Kriss et al.|[T983h and 
subtracted from the spectrum extracted from that annulus. In this 
way the deprojection routine moves inwards subtracting the con- 
tribution of projected spectra from each successive annulus to pro- 
duce a set of deprojected spectra. 

The deprojected spectra were also analysed in the energy 
range 0.5 — 7.0 keV and grouped with a minimum of 30 counts 
per spectral bin. Appropriate response and ancillary response files 
were generated as before. The deprojected spectra were fitted in 
XSPEC with an absorbed MEKAL model and the parameters set as 
previously described for the projected spectra. 

4.3 NW sector: upstream shock 

Figs. [8] and |9] show the projected and deprojected radial profiles 
for the NW sector of Abell 2146. The projected temperature is ap- 
proximately constant at ~ 6.5 keV through the length of the ram 
pressure stripped trail of gas from the cool-core. Beyond a radius 



of 50 arcsec (200 kpc) from the central AGN, marking the approx- 
imate end of the subcluster tail, the temperature increases steadily 
to a peak of 13'^2 keV at 100 arcsec (400 kpc). Then at the surface 
brightness edge the temperature plummets down to 4.5^07 keV. 
The metallicity parameter was not well constrained by the spec- 
tral fits and the increase seen at 100 arcsec is not significant. The 
sharp drop in the temperature at this radius is matched by a drop 
in the electron density shown by the deprojected profiles in Fig. |9] 
confirming that this is a shock. 

The deprojection routine DSDEPROJ assumes spherical sym- 
metry, which may not be reasonable given the cluster's highly elon- 
gated morphology in this sector. However, the steep gradient of the 
surface brightness profile outside 40 arcsec radius (Fig.O reduced 
the significance of the projected outer layers. The deprojected tem- 
peratures were therefore consistent with the projected temperatures 
to within the lC7 errors (Figs. [6] and (Sj. We used the deprojected 
results to analyse the shock at 120 arcsec where the surface bright- 
ness contours tend toward circular at the shock edge. The elongated 
morphology and substructure of the subcluster tail, coupled with 
the shallow surface brightness gradient, made the deprojection of 
the surface brightness edge at 50 arcsec much more difficult. For a 
good constraint on the temperature, a much greater number of X- 
ray counts is needed to facilitate the use of smaller radial bins in a 
narrower NW sector. The metallicity was poorly constrained in the 
deprojected annuli therefore this parameter was fixed to the average 
of 0.4 Zp determined from the spectral fits to the projected annuli. 

There are two sharp drops in the deprojected electron density 
profile (Fig. |9ll which correspond to the outer surface brightness 
edge at 120 arcsec seen in Fig.|5](right) and the subcluster tail end 
at 50 arcsec. At a radius of 120 arcsec, the density decreases by a 
factor of 2. 83 ±0.08, which coincides with a temperature drop from 
16^3 keV down to 4.6^07 keV confirming that this outer edge is a 
shock. By combining the deprojected temperature T and electron 
density ng, we calculated the electron pressure P,, = kgneT (Fig.|9] 
bottom). As expected, there is a large decrease in pressure by factor 
of at the shock front. 

The nature of the edge at the subcluster tail end (50 arcsec) 
was more difficult to determine. The radial temperature bin outside 
the edge (50 — 70 arcsec) is likely to contain some hotter gas from 
the outer shock edge in projection, biasing the projected tempera- 
ture value high. The deprojected temperature in this radial bin is 
poorly constrained so it is unclear whether the subcluster tail ends 
in a shock or a cold front. The density decreases by a factor of 
1.91 ±0.07 producing a drop in the electron pressure by a factor of 
2. 1+0 g. The decrease in pressure across the surface brightness edge 
was not very significant primarily because of the poor constraints 
on the gas temp erature. 

Following Landau & Lifshit3 ( Il959l) , we applied the Rankine- 
Hugoniot jump conditions across the shock at 120 arcsec radius to 
calculate the Mach number, M = v/c^, where v is the velocity of the 
pre-shock gas and c, is the velocity of sound in that gas. The Mach 
number can be calculated independently from the density jump, 

/ 2 P2 /. \ 1/2 

"H r±i-4Vi) ) 

or temperature jump, 

mJ(Z±14MV'^ (2) 
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Figure 8. NW sector projected radial temperature (upper) and metallicity 
(lower) profiles centred on the AGN. The black circles show the radial tem- 
perature profiles with the metallicity left as a free parameter and the red di- 
amonds show the results when the metallicity was fixed to 0.4 Zq (slightly 
offset in the x-direction for clarity). 



where T\, pi and T2, P2 denote the temperature and density be- 
fore/upstream and after/downstream of the shock respectively. Here 
we assume the adiabatic index for a monatomic gas, y = 5/3. This 
may not be applicable if, for example, a significant amount of en- 
ergy is lost in the acceleration of particles at the shock front. 

At the NW shock, the density drops by a factor P2/P1 = 
2.83 ±0.08 which, from equation [T] gives a Mach number M = 
2.7 ±0.1. Using equation |2l the observed temperature drop of 

These two independent cal- 



0.5 



T2/T1 = 3.4+i:^ gives M = 2.T 
culations of the Mach number agree within the IcT errors. 

However, the sharp drop in surface brightness produced by the 
shock edge is superimposed on the underlying decay in the clus- 
ter surface brightness with radius. The relatively large radial bins 
used in this analysis to ensure a reliable calculation of temperature 
may therefore overestimate the density drop at the shock. If free- 
free emission dominates, the X-ray emissivity depends strongly on 
gas density and only weakly on temperature, e oc p^j" 1/2 permit- 
ting narrower radial bins for calculating the density. Since the ICM 
temperature in Abell 2146 only drops below 2 keV in the cool sub- 
cluster core this is a reasonable assertion. We therefore deprojected 
the X-ray surface brightness profile, using small corrections for the 
temperature variation, and derived the electron density in narrower 
radial bins to more accurately calculate the density jumps across 
the shock edges. 

The surface brightness profile was deprojected by assuming 
spherical symmetry and using the straightforward 'onion-peeling' 
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Figure 9. NW sector deprojected radial temperature (upper), electron den- 
sity (centi'e) and electron pressure (lower) profiles centred on the AGN. 
The metallicity was fixed to 0.4Z, .j. The projected temperature profile, with 
metallicity also fixed to O.AZ ;,, is shown overlaid on the deprojected pro- 
file (slightly offset in the x-direction for claiity). The projected electron 
pressure was calculated using the projected temperatures and deprojected 
electron densities. The dashed line marks the tail end of the subcluster and 
the solid line marks the upstream shock. 



method first described in lFabianetal.l ( ll98(]|) . The background was 
subtracted using the region at large radii from the cluster as before 
(Fig.|4j. The radial bins are slightly larger for the deprojected pro- 
file to ensure a similar number of counts at each radius. However, 
there were few data points above the background level outside the 
NW shock edge to constrain the upstream gas properties. Although 
there are ROSAT observations of Abell 2146, in the PSPC expo- 
sures of the field it is more than 40 arcmin off axis and the HRI 
exposure of the cluster on axis does not detect the cluster outskirts 
with any significance. It was therefore not possible to verify the 
analysis of the outer cluster gas layers with ROSAT data. 

The deprojected electron density was calculated from the 
XSPEC normalization for a PHABS(MEKAL) model, where the to- 
tal model flux was set equal to the deprojected surface brightness 
in each radial bin. The temperature parameter was set to the corre- 
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Figure 10. Upper: Surface brightness profile in the energy range 0.3- 
7.0 keV for the NW shock. Lower: deprojected electron density profile. The 
shock edge is marked with a solid line. 



spending deprojected value determined from spectral fitting (Fig.|9] 
upper). 

The projected surface brightness and deprojected electron 
density across the NW shock edge are shown in Fig. [TO] Fig. [TO] 
(lower panel) shows a drop in the gas density of pi/pi = 2.0 ±0.4 
across the shock edge. From equation[T] this gives a Mach number 
M = 1 .7 ± 0.3, lower than the spectral fitting result calculated from 
the density jump of M = 2.7 ± 0. 1 . If we repeat the spectral depro- 
jection, but with each radial bin split into two and the temperature 
parameters tied together in each pair, the Mach number calculated 
from the density also drops to M = 2. I^q j . 

We therefore conclude that the Mach number of the NW 
shock is closer to the lower value of M = 1.7 ±0.3 calculated 
from the higher resolution spatial binning which more closely 
probes the shock. The sound speed in front of the shock is Cj = 
{ykBTi/mH^iYl^ = llOO+gy" km s^' , where the gas temperature 

1.0 



in front of the shock T\ = 4.6_g -j 



keV and the mean molecular 
weight of the medium jl — 0.6. Therefore, for a Mach number of 
M = 1 .7 ± 0.3 the shock velocity is v = Mcs = 1900 ± 400 km s" ' . 



4.4 SE sector: cold front & bow shock 

The projected and deprojected radial profiles in temperature, den- 
sity and electron pressure for the SE sector in front of the subcluster 
core are shown in Fig. [TT] The assumption of spherical symmetry 
applied by DSDEPROJ was considered to be reasonable as the clus- 
ter appears approximately circular on the sky within this sector and 



the radial bins traced the surface brightness contours (Fig.|7]l. The 
metallicity was poorly constrained in the spectral fits and so was 
fixed to an average value of 0.4 Zq in this sector. 

The subcluster core, analysed by the innermost radial bin, con- 
tains the lowest temperature, 2.4 ±0.1 keV, and highest density 
gas, = 0.0130 ±0.0001 cm^-^ , in the cluster. The radiative cool- 
ing time of the gas J^ool was derived from the temperature T and 
density iie 



^cool ' 



5 nkgT 
2 nenHA(T) 



(3) 



where A(r) is the cooling function, n is the total number density of 
gas particles and is the number density of hydrogen. In the sub- 
cluster core the radiative gas cooling time drops to 3.5 ± 0.2 Gyr. 
Therefore it is likely to be a cool core remnant that is being ram 
pressure stripped in the merger. 

The temperature of the gas increases rapidly in front of the 
subcluster core by a factor of 1 -S^q ^. This jump coincides with the 
surface brightness edge at 18 arcsec visible in Fig.[5]and a drop in 
the electron density by a factor of 3.3^q[. The electron pressure 
is therefore approximately constant across this surface brightness 
edge. This edge is a cold front or contact discontinuity, similar to 
that found at the leading edge of the bullet in the Bullet cluster 
( Markevitch et al. 2002 ) and in other cl usters such as Abell 2 142 
i Markevitch etal iboool) and Abell 3667 iVikhhnin et"aDl200lh . 

The gas properties at the outer edge at 55 arcsec radius were 
more difficult to extract. The sharp drop in the cluster surface 
brightness necessitated a large outermost radial bin, which there- 
fore did not allow constraints on the properties close to the edge. 
There is a very sharp drop in the electron pressure by a factor of 
18+^' at 55 arcsec corresponding to a decrease by factors of 3_j 
and 6.0^0 2 in temperature and density respectively. However, the 
Rankine-Hugoniot jump conditions (eq. \Q give a maximum pos- 
sible density contrast across adiabatic shocks in a monatomic gas 
of a factor of four. This limit arises because as the strength of the 
shock is increased to higher Mach numbers the high thermal pres- 
sure building behind the shock limits the compression of the post- 
shock gas. The large radial bins (40 ±11 arcsec and 93 ± 42 arcsec) 
used for this analysis do not closely probe the gas properties across 
the shock edge. The underlying decline in the cluster density pro- 
file, which is significant in such a large radial region, is included in 
the calculation of the density change at the shock. This produces an 
overestimate of the density jump. 

We have therefore also calculated the density jump by de- 
projecting the surface brightness profile across the shock. Fig. [12] 
shows the drop in surface brightness across the shock edge at 
55 arcsec radius corresponds to a drop in the gas density by a much 
lower factor of pa/pi = 2.4±0.7. Using equation[T] the Mach num- 
ber of the SE bow shock is therefore M = 2.2 ± 0.8. For a preshock 
temperature of 4^^ keV, which gives a sound speed of cs = 



1000+^i]j] km s- 



the shock velocity is v = 2200 



-1000 

-yoo 



km s 



The subcluster is at a projected distance of roughly 380 kpc from 
the centre of the main cluster and, if we assume that the subc luster 
moves at the shock velocity (but see lSpringel & Farradl2007l) . this 
implies that the subcluster passed through the centre of the main 
cluster approximately 0.1—0.3 Gyr ago. 
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Figure 11. SE sector projected and deprojected radial temperature (upper), 
electron density (centre) and electron pressure (lower) profiles centred on 
the AGN. The metallicity was fixed to OAZ,: . The projected electron pres- 
sure was calculated using the projected temperatures and deprojected elec- 
tron densities. The projected data points are slightly offset in the x-direction 
for clarity. 
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Figure 12. Upper: Surface brightness profile in the energy range 0.3- 
7.0 keV for the SE shock. Centre: deprojected surface brightness profile 
Lower: deprojected electron density profile with bestfit values on either side 
of the shock edge. The dotted line marks a possible density enhancement in 
the stagnation region in front of the subcluster 



5 DISCUSSION 

By analysing the temperature and density of the gas on either side 
of the surface brightness edges, we have determined that the merg- 
ing cluster Abell 2146 contains two shock fronts: a bow shock in 
front of the subcluster with M = 2.2 ±0.8 and a slower upstream 
shock behind the main cluster with M = \.l ±0.3. We estimate 
that the subcluster passed through the centre of the main cluster ap- 
proximately 0.1 — 0.3 Gyr ago and is being ram pressure stripped 
of its material. The sharpness of the observed surface brightness 
edges (Fig.|5} suggests that Abell 2146 is oriented close to the plane 
of the sky. However, a measurement of the line of sight velocities 
for the subcluster and the main cluster galaxies using optical spec- 
troscopy, combined with the Mach numbers, will give a quantitative 
constraint on the inclination of the merger axis. 

Although Abell 2146 is a smaller system with a much lower 
global tempera ture, it appears remarkably similar in structu re to the 
Bullet cluster l lMarkevitch et alj|2002l : lMarkevitcril2006h . Shock 
fronts are so rarely detected because the merger must be observed 



before the shock has moved to the outer, low surface brightness re- 
gions of the cluster and with a merger axis close to the plane of the 
sky so that projection effects do not conceal the surface brightness 
edges. It is therefore unsurprising that Abell 2146 appears similar 
to the Bullet cluster and is observed at a comparable time since core 
passage of the subcluster (both around ~ 0.2 Gyr). 

By choosing the brightest clusters to observe we will also 
preferentially select merging clusters in this short window of core 
passage. Merging clusters undergo a dramatic increase in lumi- 
nosity and temp erature as the cores col l ide, pushing them up the 
Lr — T relation jRicker & Sarazirj|200ll : iRitchie & Thomasll2002l: 
iMarkevitcriboogT This boost is caused by the compression of the 
cluster cores during the merger and can be up to a factor of 10 
for the luminosity in a head-on collision between two equal mass 
clusters. Statistically more likely collisions between clusters with 
higher mass ratios of 3:1 and 8:1 and low impact parameters (less 
than a few core radii) will still produce an inc rease in luminos- 
ity by a factor of ~ 4 and ~ 2, respectively l lRicker & SarazinI 
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I2OOII ; iRitchie & Thomasll2002l : IPoole et al.ll2007l) . However, this 
strong variation is relatively short-lived and disappears after ap- 
proximately one sound crossing time. For a typical cool-core clus- 
ter, with a core radius of 100 kpc and temperature 3 keV, the sound 
crossing time is only 0. 1 Gyr. 



5.1 Mass ratio of the merging clusters 

The key morphological difference between Abell 2146 and the Bul- 
let cluster is the clear detection of an upstream shock. An upstream 
shock is generated as the gravitational pot ential minimum fluc- 
tuates rapidly duri ng core passage (see eg. iRoettiger et alj[l997l : 
iGomez et all2002l) . The gravitational potential reaches an extreme 
minimum as the two cluster cores coalesce. This causes a signifi- 
cant amount of the outer cluster gas to flow inwards. The subcluster 
then exits the main core and the gravitational potential rapidly re- 
turns to its premerger level, which expels much of the newly arrived 
gas. The gas that is expelled in the upstream direction then collides 
with the residual infall from the subcluster and forms a shock prop- 
agating in the opposite direction to the subcluster. 

For a higher mass ratio merger, such as the Bullet cl uster 
(10:h lClowe et alj2004lBradac et alj2006l : lciowe et alj2006l) . the 
perturbation to the gravitational potential caused by the merger is 
smaller relative to the total potential and shorter in duration. The 
outer gas layers have less time to respond to the change in gravita- 
tional potential, which reduces the infall and strength of the sub- 
sequent outflow, producing only a weak upstream shock. In ad- 
dition, there is relatively little residual infall from the wake of 
the small subcluster. Simulations of cluster mergers over a range 
of mass ratios have shown that the more dramatic variation in 
the grav itational potential during lower mass ratio mergers (less 



than 8:1; Roettiger et al prod uces a stronger upstream sho ck 

( IRoettiger et al.lll997l:lRicker & Sarazii]|200ll : IPoole et al.ll2006h . 

The comparable strength of the bow and upstream shocks in 
Abell 2146 suggests that this mass ratio could be lower than 4:1. 
Closer to an equal mass ratio is unlikely because the predicted 
peak Mach number at core passage will then drop below the ob- 
served M ^ 2. We therefore argue that the subcluster and main 
cluster components of Abell 2146 are likely to have a mass ratio 
of 3 or 4:1. New Subaru Suprime-Cam observations of this cluster 
(PI Gandhi) will allow a more quantitative analysis of the mass. 



5.2 Subcluster structure 

As the subcluster passed through the core of the main cluster, its 
leading edge was compressed and swept back as the strong gravi- 
tational potential stripped away the gas. The local gas was pushed 
aside during core passage producing an elongation of the gas distri- 
bution perpendicular to the merger axis. For high mass ratio merg- 
ers, a ring of compressed gas forms when the small subcluster 
passe s through the primary core ( IRoettiger et alJlT997l : |Poole et al.l 
l200d) . This can be clearly seen as a bar-like structure in the obser- 
vations of the Bullet cluster ( Markevitch 2006). Figs.[3]and|6]show 
an extended spur of cool gas at 5 — 7 keV to the SW of the subclus- 
ter tail. This is likely to be stripped material from the subcluster and 
swept up gas from the main cluster. However, there does not appear 
to be a symmetric feature to the NE of the subcluster tail, although 
there is a suggestion of a smaller spur to the NE in Fig.jSj Simula- 
tions suggest that elongation perpendicular to the merger axis may 
be less evide nt in lower mass ratio mergers where t he cores are 
similar sizes iRoettiger erZlll997l ; [Poole et al.ll2006l) . A slightly 



off axis merger will also produce an asymm etrical distribution of 
swept up and stripped gas JPoole et al.ll2006h . 

The bow shock visible in front of the subcluster core formed 
when the subcluster's infall velocity exceeded the sound speed in 
the ambient cluster gas, cs = IOOOI^q" kms^' . In principle, the 
Mach cone angle should be directly related to the Mach num- 
ber of the shock. However, as has been found for Bu llet cluster 
dMarkevitch et al.ll2002l : iMarkevitch & VikhUniijboOTi) . the sub- 
cluster is shrinking over time and decelerating in the gravitational 
potential of the main cluster. Therefore, in practice the subcluster 
can not be approximated as a solid body moving at constant veloc- 
ity. ISpringel & Farraj i2007h also found in their simulations of the 
Bullet cluster that the opening angle of the Mach cone was much 
wider than expected and cannot be easily used to independently de- 
termine the shock strength. For Abell 2146, we can only estimate 
the Mach cone angle from the narrow section of the bow shock that 
is easily discernable but it is clear that it is much broader than the 
tj) — 30° expected for a M ~ 2 shock (Fig.[T). 

We might also expect that the stand-off distance between the 
bow shock and the subcluster is relate d to the Mach n umber, ac- 
cording to the approx imate relation of iMoeckel jl949lPI (see also 
IVikhlinin et"aDl200lh . This analysis requires measurements of the 
geometry of the projectile, which was particularly difficult to de- 
termine for the subcluster gas cloud. We estimated the half-width 
of the subcluster core to be ^ 35 kpc, which for M = 2.2 ± 0.8 
gives a predicted stand-off distance between the shock and the sub- 
cluster of d ~ 20 — 40 kpc. The observed distance is approximately 
140 kpc. A similar result can be obtained for the Bullet cluster 
The assumptions of a uniform preshock medium, constant veloci- 
ties and a solid projectile are clearly not applicable for these merg- 
ing subclusters. 

Behind the subcluster, the ambient cluster gas that was pushed 
aside during its passage will fall back and produce tail shocks (eg. 
IRoettiger et alJl997l : |Poole et al.l200^ . Fig.|9|shows a drop in den- 
sity by a factor of 1.91 ±0.07 at the approximate position of the tail 
end (--^ 50 arcsec NW from the AGN). If this feature were a shock 
we would expect a corresponding jump in the temperature by a 
factor of 1.66 ±0.07. Despite the uncertainties in the temperature 
values (section |43] | such a large decrease appears unlikely. Alter- 
natively, this feature could be a cold front separating the cooler, 
ram pressure stripped material from the subcluster which is slowly 
falling back into the hotter main cluster gas. Although there is a de- 
crease in the electron pressure across the front, the large uncertain- 
ties in the temperature value reduce the significance of this drop. 
A longer Chandra observation providing a more accurate measure- 
ment of the temperature could resolve this issue and allow con- 
straints on the ICM transport processes across the shocks a nd cold 
fronts (eg. lMarkevit"ci3l2006l : IMarkevitch & Vikhlininll2007l) . 

The subcluster will also generate strong turbulence in its wake. 
Figs. [T] and |2] show complex structures in the ram pressure stripped 
material of the subcluster tail. The temperature map also suggests 
there could be substructures and smaller shocks inside the tail. The 
ram pressure stripped tail is considerably warmer (5 — 8 keV) than 
the subcluster core where it originated (1.9 ±0.1 keV). This level 
of heating appears to be consistent with the results for the Bullet 
cluster where the subcluster core gas at ~ 2 keV is stripped and 
forms a tail of gas heated to 7 — 10 keV iMarkevitch et alj|2002l : 
iMillion & AUenl I2OO9I) . Along the surface of the bullet, Kelvin- 
Hehnholtz and Rayleigh-Taylor instabilities are expected to de- 



^ Available at |http://naca.central.cranfi eld.ac.uk/reports/1949/naca-tn- 1921.pdf | 
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Figure 13. Emission measure per unit area map (units cm ^ arcsec ^ ). 
The emission measure is the XSPEC normalization of the MEKAL spectrum 
K = EII{A X \0^*nDj^(\+zf ), where EI is the emission integral EI = 



velo p and break up the subc luster, as has been observed in Abell 
520 jMarkevitch e"taDl2005h . The timescale for the destruction of 
the subcluster will also likely be influenced by magnetic fields 
which sta bilize against these instabilities (eg. Ijones et al] 1 19961 : 
IVikhlinin et a l. 2001). Fig. [13] shows an emission measure map of 
the cluster produced with the Contour binning algorithm (section 
|4.1| l but with a lower signal-to-noise ratio of 15 giving better spa- 
tial resolution at the expense of larger errors (~ 15%). The absorp- 
tion parameter was fixed to the Galactic value and the abundance 
was fixed to O.4Z0 for the spectral fitting. The finer spatial binning 
in Fig. [T3]hints at possible substructure in the subcluster tail but a 
deeper Chandra observation is required to significantly detect any 
structure. 



5.3 Location of the brightest cluster galaxy 

In a merger event, the galaxies are effectively collisionless particles 
and therefore, along with the dark matter component, lead the X- 
ray gas on exiting the main cluster core. The X-ray gas is slowed 
as a result of the ram pressure from the intera cting cluster cores . 
This is clearly observed in the Bullet cluster (eg. lClowe et al .l2006h . 
However, Fig. |2] shows that while there are many galaxies leading 
the X-ray subcluster core in Abell 2146, the BCG is located imme- 
diately behind the X-ray peak. 

Observations of the merging cluster Abell 168 also found 
that the subcluster cD galaxy lagged behind the cool gas peak 
jHallman & MarkevitchI 120041) . As the subcluster passes through 
the core of the primary, ram pressure pushes back the gas from 
the gravitational potential. Then as the subcluster enters the outer, 
less dense cluster gas, the ram pressure drops rapidly and the 
cool gas core rebounds and oversho ots the subcluster dark matter 
peak in a 'ram pressure slingshot' faallman & MarkevitclJ[2004 



iMathis et alj2005l:lAscasibar & Markevitcril2006l) . The cD galaxy, 
which should trace the local gravitational potential, will then ap- 
pear behind the subcluster' s cool core. However, this represents a 
late stage of the merger where the subcluster has reached its apoc- 
entre. For A bell 2146, this wo uld occur ~ 1 Gyr (based on a mass 
ratio of 3: 1, iPoole et alj|200^ after the subcluster passes through 
the main cluster core and is therefore inconsistent with the esti- 
mated 0.1 — 0.3 Gyr for the age of the merger. Although it was 
difficult to determine the location of the main cluster core from the 
existing observations, the prominence of the two shocks and the 
undisrupted subcluster core provides strong evidence against this 
being a late stage merger. Unambiguous detections of shock fronts 
are so rare because they can only be seen in the early stages of a 
merger, before the shock front h as propagated to the low surface 
brightness outskirts of the cluster jMarkevitch & Vikhlininll2007i) . 

The HST ob servation of Abell 2146 (contours shown in Fig.|2j 
ISand et aill2005l) shows that the BCG has a large halo of emission, 
which is extended along the merger axis and clearly traces the X- 
ray emission in the subcluster. The large diffuse envelope of the 
BCG may therefore be interacting with the surrounding X-ray gas. 
The BCG does not obviously appear to be disrupted by the merger 
event so it remains unclear why it is located in the wake of the 
subcluster cool core. 

In an imaging survey with the Spitzer Space Telescope of 62 
BCGs with optic al line emission loca ted in the cores of X-ray lu- 
minous clusters, lOuillen et alj ( |2008|) found that Abell 2146 has 
a high IR luminosity and the secon d highest rate of st ar forma- 
tion in the sample at 192M0 yr"' . lO'Dea et al] j2008|) found a 
correlation between the mass deposition rates estimated from X- 
ray observations of the sample and the IR star formation rates, 
where the star formation rate is ~ 1/10 of the mass deposition rate. 
This suggests that the cooling ICM is the sour c e of the gas that 
is forming stars (see also Johnst one et alJI 19871: iMcNamara et al 



2004lHicks & Mushot zkv 2005: IRaffertY et al.l2006l : ISalome et ai 



2006h . However, in Abell 2146, the high star formation rate of 



192Mq yr~ is ~ 5 times greater than our estimated mass depo- 
sition rate of 40 ± lOM© yr^^ . This high rate of star formation in 
the BCG could have been triggered by the cluster merger and this 
will be explored in a future paper. 



5.4 Limits on diffuse radio emission 

The connection between cluster mergers and the presence of dif- 
fuse, steep spectrum radio sources has been extensively inve sti- 
gated (Tribblel 993l : lRoettiger eTakll 19991 : ISninetti et al.H2009l) so 
we have searched the radio survey data around Abell 2146 to de- 
termine a limit to any diffuse emission. Unfort unately, Abell 2146 
is just below the X-ray luminosity limit set bv lVenturietal]|2007l 
for their 610 MHz GMRT survey so no targetted radio imaging 
at frequencies below 5 GHz exists. However, if we take the ob- 
served upper bound for radio halo power from .Brunetti et al ] |2009l 
for a cluster with the X-ray luminosity of Abell 2146 (Pi.4ghz ~ 
1q24.0 Hz^l ) we would expect a halo or relic that is weaker than 
~6mJy at 1.4 GHz. 

Inspecting the VLSS 74 MHz, WENSS 327 MHz and NVSS 
1.4 GHz survey data and archival VLA databases, we find two un- 
resolved sources detected in WENSS, NVSS and 5 GHz VLA im- 
ages. One is coincident with the BCG and is detected at 15.3 ± 
0.6 mJy at 1.4 GHz with a spectral index of -0.43. The other is 
associated with a probable cluster member that is 40.6 ±1.3 mJy 
at 1.4 GHz with a spectral index of -0.48. These two sources are 
blended in WENSS catalog but the total flux density is consis- 
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tent with the sum of the two point source components. Inspect- 
ing the noise in these maps we calculate 3(7 flux density limits of 
< 660 mJy, < 21 mJy and < 2.5 mjy at 74, 327 and 1420 MHz 
respectively for any diffuse emission on scales of 1 — 2.5 arcmin. 
Therefore, a radio halo a factor of around two below that expected 
from other comparable systems is consistent with the current obser- 
vations irrespective of the spectral index. A deep eVLA or GMRT 
observation is required to improve on this limit. 



6 CONCLUSIONS 

The Chandra observation of Abell 2146 has revealed a merg- 
ing system where a ram pressure stripped subcluster has recently 
passed through and disrupted the primary cluster core. From 
the X-ray temperature and surface brightness maps, we found a 
bow shock propagating in front of the cool 2 — 3 keV subclus- 
ter and calculated a Mach number M — 2.2 ± 0.8 from the den- 
sity jump across the shock. The subcluster velocity is therefore 
V = 2200+500" ^® estimated that the subcluster passed 

through the main cluster core 0.1 — 0.3 Gyr ago. In addition, there 
is a factor of IO+2 drop in the electron pressure in the outskirts 
of the main cluster indicating the presence ofaM = 1.7±0.3 up- 
stream shock. There are potentially further shocks and a cold front 
in the gas tail behind the subcluster but these features could not be 
confirmed with the existing data. 

Although Abell 2146 is a smaller and cooler system than the 
Bullet cluster, it appears similar in structure and to be at a compa- 
rable merger epoch of ~ 0.2 Gyr since core passage. This can be 
understood as a selection effect: shock fronts can only be detected 
at this early stage in the merger evolution before they have prop- 
agated to the outer, low surface brightness regions of the cluster. 
The merger axis must also be close to the plane of the sky so that 
projection effects do not conceal the surface brightness edges. 

Based on the measured shock Mach numbers of M 2 and the 
strength of the upstream shock, we estimate a mass ratio between 
the two merging clusters of around 3 or 4:1. Forthcoming Subaru 
observations of Abell 2146 will allow a more quantitative analysis 
of the mass distribution between the two clusters. We compared the 
Chandra observation with an archival HST observation and found 
that while there is a group of galaxies located in front of the X-ray 
subcluster core, the brightest cluster galaxy is located immediately 
behind the X-ray peak. A future weak lensing analysis coupled with 
galaxy velocities along the line of sight from optical spectroscopy 
could help to explain the galaxy cluster dynamics. 
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